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Abstract 

In this work we investigate the background dynamics when dark energy is coupled to dark matter with 
a suitable interaction in the universe described by Loop quantum cosmology. Dark energy in the form of 
Generalised Cosmic Chaplygin gas is considered. A suitable interaction between dark energy and dark 
matter is taken into account in order to at least alleviate (if not solve) the cosmic coincidence problem. 
The dynamical system of equations is solved numerically and a stable scaling solution is obtained. A 
significant attempt towards the solution of the cosmic coincidence problem is taken. The statefinder 
parameters are also calculated to classify the dark energy model. Graphs and phase diagrams are drawn 
to study the variations of these parameters. It is seen that the background dynamics of Generalised 
Cosmic Chaplygin gas is completely consistent with the notion of an accelerated expansion in the late 
universe. From the graphs, generalised cosmic Chaplygin gas is identified as a dark fluid with a lesser 
negative pressure compared to Modified Chaplygin gas, thus supporting a 'No Big Rip' cosmology. It 
has also been shown that in this model the universe follows the power law form of expansion around the 
critical point, which is consistent with the known results. Future singularities that may be formed in 
this model as an ultimate fate of the universe has been studied in detail. It was found that the model is 
completely free from any types of future singularities. 

1 Introduction 

At the turn of the last century observations from la Supernova confirmed that our universe is suffering from an 
accelerated expansion [TJ [2] . In the quest of finding a suitable model for the expanding universe, Cosmologists started 
to investigate the root cause that is triggering this expansion. Fundamentally, we were to modify Einstein's equation 
either by modifying the left hand side, i.e., modifying the idea of Einstein gravity or to modify the right hand side 
which immediately speculate the nature of the matter inside the universe. If our Universe is filled up by some invisible 
fluid having a large negative pressure then it violates the strong energy condition i.e. p + 3p < 0. Due to its invisible 
nature this energy component is aptly termed as dark energy (DE) [3]. 

With the introduction of DE, the search began for different candidates that can effectively play the role of DE. 
DE represented by a scalar field |f| is often called quintessence. Not only scalar field but also there are other Dark fluid 
models like Chaplygin gas which plays the role of DE very efficiently. The earliest form of this was known as pure 
Chaplygin gas [H |5] . Extensive research saw pure Chaplygin gas modify into generalized Chaplygin gas [H 13 El [10] . 
Finally modified Chaplygin gas (MCG) came into existence [HJH2]- Dynamics of MCG in Loop quantum cosmology 
was studied by Jamil et al [13]. Dynamics of MCG in Braneworld was studied by Rudra et al [14j . 

Currently, we live in a special epoch where the densities of DE and dark matter(DM) are comparable, although 
they have evolved independently from different mass scales. This is known as the famous cosmic coincidence problem, 
till date several attempts have been made to find a solution to this problem [15 l fl6 l [17 ] rTH I ri9 t [20 ]. [21 ^ [22 | [23] , f24 \ [25 ] . 
A suitable interaction between DE and DM is the best possible tool to find an effective solution to this problem. It 
is obvious that there has been a transition from a matter dominated universe to dark energy dominated universe, by 
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exchange of energy at an appropriate rate. Now the expansion history of the universe as determined by the supernovae 
and CMB data [2TJ [22] bounds us to fix the decay rate in such a way that it is proportional to the present day Hubble 
parameter. Keeping the above fact in mind cosmologists all over the world have studied and proposed a variety of 
interacting DE models [25J 123 GH [Ml 1301 EU 132] ■ 

In 2003, P. F. Gonzalez-Diaz [33] introduced the generalized cosmic Chaplygin gas (GCCG) model. The speciality 
of the model being that it can be made stable and free from unphysical behaviours even when the vacuum fluid 
satisfies the phantom energy condition. In the previous studies related to DE corresponding to phantom era Big-Rip 
is essential, as the time gradient of scale-factor blows to infinity in finite time. For the first time P. F. Gonzalez-Diaz 
including the GCCG model showed that Big Rip, i.e., singularity at a finite time is totally out of question. Hence in 
such models there is no requirement for evaporation of black hole to zero mass. The Equation of state (EOS) of this 
model is 

(1) 



P = -P~ U 



c + | p (i+«)_cV 



where C = — 1, with A being a constant that can take on both positive and negative values, and — C > u> > 0, C 
being a positive definite constant, which can take on values larger than unity. GCCG can explain the evolution of the 
universe starting from the dust era to AC'DM, radiation era, matter dominated quintessence and lastly phantom era 

As we have stated earlier, modifying the right hand side of Einstein's equation was not the only way to explain the 
increase in the rate of expansion of the universe. The gravity part on the left hand side of Einstein's equation can also 
be modified to demonstrate the present day universe. In this perspective Loop quantum gravity (LQG) was introduced 
and Loop quantum cosmology (LQC) was developed. In this work we consider LQG as the modified gravity theory. 
Our main aim of this work is to examine the nature of the different physical parameters for the universe around the 
stable critical points in LQC model in presence of GCCG. Impact of any future singularity caused by the DE in LQC 
model will be studied as well. 

This paper is organized as follows: Section 2 comprises of the basic equations of the LQC model. Section3 
deals with the dynamical system analysis of GCCG in LQC model. In Section 4 we present a complete graphical 
representation of phase plane analysis. In Section 5 we study the future singularities, followed by some concluding 
remarks in section 6. 



2 The Model: Loop Quantum Gravity 

In recent years, loop Quantum Gravity (LQG) has evolved as an outstanding effort to describe the quantum effect 
of our universe [35l [36]. LQG is a theory trying to quantize the gravity with a non-perturbative and background 
independent method. The theory and principles of LQG when applied in the cosmological framework creates a new 
theoretical framework of Loop Quantum Cosmology (LQC) 37, 38, 39 . In this theory, classical space-time continuum 
is replaced by a discrete quantum geometry. The effect of LQG can be described by the modification of Friedmann 

2 

equation by adding a term quadratic in density. In LQC, the non-perturbative effects lead to correction term ^ to 
the standard Friedmann equation. With the inclusion of this term, the universe bounces quantum mechanically as the 
matter energy density reaches the level of pi (order of Plank density). 



2.1 Basic equations of Loop quantum gravity model 

Recently the model of DE has been explored in the framework of LQC. The cosmological evolution in LQC has been 
widely studied for quintessence and phantom DE models |40j . The modified Friedmann equation for Loop Quantum 
Cosmology is given by [2"9"I |4"TI |4"2"] 

h 2 = Pt( 1 -8Z] ( 2) 



Pi . 

Here p\ = "v/37r 2 7 3 G 2 ?i is the critical loop quantum density and 7 is the dimensionless Barbero-Immirzi parameter. 
Pt — Pm + Pe represents the total cosmic energy density, which is a sum of energy density of DM (p m ) and the energy 
density of DE (pe)- 

As in the present problem the interaction between DE in the form of GCCG and pressureless DM has been taken 
into account. For interacting GCCG and DM the energy balance equation will be 

Pgccg + 3H(1 + w gccg ) p gccg = -Q, for GCCG (3) 
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and 

p m + 3Hp m = Q, for the DM interacting with GCCG. (4) 

where Q = 3bHp is the interaction term, b is the coupling parameter (or transfer strength) and pr = p = Pgccg + Pm 
is the total cosmic energy density which satisfies the energy conservation equation p + 3H (p + p) = [4"31 [T3] . 

As we have lack of information about the fact, how DE and DM interact so we are not able to estimate the 
interaction term from the first principles. However, the negativity of Q immediately implies the possibility of having 
negative DE in the early universe which is overruled by the necessity of the second law of thermodynamics to be held 
[44] . Hence Q must be positive and small. From the observational data of 182 Gold type la supernova samples, CMB 
data from the three year WMAP survey and the baryonic acoustic oscillations from the Sloan Digital Sky Survey, it 
is estimated that the coupling parameter between DM and DE must be a small positive value (of the order of unity), 
which satisfies the requirement for solving the cosmic coincidence problem and the second law of thermodynamics 
[15] . Due to the underlying interaction, the beginning of the accelerated expansion is shifted to higher redshifts. 
Consequently using the Friedmann equation ([2]) and the conservation equation, we obtain the modified Raychaudhuri 
equation 



H = -\{ P + P) (1-2^ 
2 V Pi 



(5) 



3 Dynamical system analysis 

In this section we plan to analyse the dynamical system. For that firstly we convert the physical parameters into some 
dimensionless form, given by 

x = hxa, u=— T , v = —; (6) 



3H 2 ' 3iJ 2 

where the present value of the scale factor, ao = 1 is assumed. Using equation ([T]) to equation ([6]) we get the parameter 
gradients as given below: 



— = 3 

dx 



l-v 
1 + v 



{V + U (1 + W gccg )} ~ b (U + V) - U (1 + W gCC g) 



and 



dv 
dx 



= 3 



b (u + v) - v + v Q — {v + u (1 + Wg CC g)} 



Where, w gcC g is the EoS parameter for GCCG determined as 



Pgc 



(u + vf 



■c-f-c+ " 2( "+r^ v 

V (u+v) ) 



gecg 



Pgccg U 2 (u + V — l) 2 p\ 

In the above calculations for mathematical simplicity we have considered a = 1. 



(7) 
(8) 

(9) 



3.1 Critical points 

The critical points of the above system are obtained by putting 4^ = = 4^. Now this system of equation does not 
yield an explicit solution. So we have to investigate for a numerical solution by putting some special values to the 
parameters which will yield a non zero positive solution of the system of equations. Considering uj = —1, we obtain 
the following critical values of the system. 



1 (\ 1 

U ^=2b 4 + 4 P - 



Q 

2V2 



and 

1/1 1 Q 



4 4 
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Vlc = I + I P "27! 
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-P 



2V2 y 
(10) 

(11) 



2V2, 

(12) 
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1 1 Q 



4 4 2^2 



(13) 



where P = y/l - 46 and Q = ./-3 
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The above two critical points correspond to the era dominated by DM and GCCG type DE and exist for 6 < j. 
For the critical points the equation of state parameter ^ of the interacting DE takes the form 



(u ic + V ic ) 



w gccg — 



-c- -c 



t4("ic+t>ic-l) 2 Pi N 

(n ic +i, ic ) 4 



(14) 



where i = 1, 2 and which holds for iti C + Uj c 7^ 1. 



3.2 Stability around critical point 



Now we check the stability of the dynamical system (eqs. (0 and ((8])) about the critical points. In order to do this, 
we linearize the governing equations about the critical points i.e., 



u = u c + Su and v = v c + Sv, 



Now if we assume / = 4 s and g = 4^, then we may obtain 



[duf]Ju+[d v f}Jv 



and 



S ' else) = [ 9 "5] c <5w+ [5„p] c 5« 



(15) 

(16) 
(17) 



where 



<9«/ = [3 {-CB 2 Af 3 ((4 + 6 - 4u) (u -l)u+(u + u 2 (4u - 11) + 6 (u 2 + u - l)) u + (1 + u (7u + 26 - 5)) v 2 

+ (Au + b - 2) u 3 + v 4 ) + -^N^ (~u 2 B 2 M 3 (2m 3 (v - 1) - 3u (u + 1) - 2u (1 + u (3 + (v - 1) u)) w 

vrp\ 

-v (v 2 - l) (1 + 2u) + u 2 (5 + v (-5 -2v(w- 1)) p 2 + u 4 B 5 (l + 2u (u - 1) + v 2 + 6 (1 + u)) p 4 A^ 



+CM 7 (2u 3 (u - 1) + « - v 3 - 3u (u + 1) + u 2 (5 + v (2v - 5))) (1 + CN U )))) 



(v + l)B 3 (CM 4 — u 2 B 2 p\) 



(18) 



d v f = 



(i + vy 



2u 2 M 4 



|-m 3 B 5 (b (v + i f + u (-1 + 2u + v (2 + 



CM 4 -ii 2 B 2 p? -CuB 3 M 4 p 2 +u 3 B 5 p 4 
p\N w - 2CM 7 (-2 + m 3 - 3v + v 3 + u (3 + v 3 ) + u 2 (-2 + v (2 + u))) (1 + CW W ) + uB 2 M 3 p\ (2u (v + l)(u + v-2) 
(1 + u (v - 1) + v) (1 + w) + C [Au 4 + b (v - 1) v (v + l) 2 + u 3 (u (3<u + 10) - 7) + u 2 {l + b (v + l) 2 + v (v + 2) (Av - 3) 



+u (-4 + 6 (« + l) 2 (2v - 1) + v (v + 1) (-5 + v (v + 2))) ) iV w ) } 



(19) 



d u g = 3[b+{v(v-l) (-u 4 B b p\N" - CM 7 (u(u-3) + v- v 2 ) (1 + CW") - u 2 B 2 M 3 p\ (-3u + u 2 + v-v 2 



-2ulj + 2voj - 2uvu - 2i; 2 u; + 2Cu (u + v- 2) N"))} 



(1 + v) B 3 p 2 {-CM 4 + u 2 B 2 p 2 ) 



(20) 



d v g = 3 



1 + 6 - u - 2v ■ 



2 (u - 1) 
(l + «) 2 



+ {M 3 7V-" (CM 4 (v (« - 1) (3u 2 - 5) + u 2 (v (v + 2) - 1) + u (1 + v (v (Av + 3) 



-6))) (1 + CN U ) - u 2 B 2 p\ (u - u 2 + 5v - 6uv + 2u 2 v ~ 5v 2 + 3uv 2 + u 2 v 2 - 3v 3 + 4uu 3 + 3w 4 + 4wcj - 2uvu - 2v 2 uj 
-4v 3 lu + 2v a lo + C(v(v-1) (3v 2 - 5) + u 2 (v (v + 2) - 1 + u (1 + v (-6 + v (4t> + 3)))) N u ))) 

1 



(1 + vf B 3 p 2 {CuM 4 - u 3 B 2 p 2 ) 



(21) 



where B = u + v — 1, M = u + v, N = —C 



3.3 Nature of cosmological parameters 

We calculate the deceleration parameter q = — 1 — (H/H 2 ), in this model as 



3 ( . Pgccg \ I jj_ 



_£_ 

' in 



q = -l + ^[l+u Jgccg ^j {T^) (22) 
Now expressing the above expression in the form of dimensionless density parameter £l gccg — Pa "p a , we get, 

3 /i_2-£-\ 
q =-! + -{! + uj gccg n gccg ) I - _ ^ J 

From the above expression it is evident that in the limit p\ — > oo, we retrieve the result for Einstein's gravity as 
follows, 



3 
2 

Now assuming = e ~ O(l) and using the transformations of equation ([5]), we get, 



QEG — ^1 + 7T (1 + Lu gcc g Q gccg ) 



Since a physically acceptable solution corresponds to both the critical points, such that (u,v) — > (ui C ,Vi C ) , i = 1,2. 
Hence from equation ([2"31 we get, 

q c = -l + l Z , where Z=(l + ^^)(^), i = 1, 2 (24) 



2 \ <; u ' !r / \ ! < / 

As special cases, we observe that for e = i, we have g = — 1 while e = 1, we have g = — oo if w gC cg > ^il - ' 
We know that the Hubble parameter varies as 

"=4 (25 > 

where we have ignored the integration constant. Integration of the above equation (|25[) yields, 

a(t)=a t& (26) 

which, as its form suggests gives a power law form of expansion. 

Sahni et al [46] introduced a pair of cosmological diagnostic pair {r, s} which is known as statefinder parameters. 
The two parameters are dimensionless and are geometrical since they are derived from the cosmic scale factor alone. 
Also this pair generalizes the well-known geometrical parameters like the Hubble parameter and the deceleration 
parameter. The statefinder parameters are given by 

r =^P' S = 3( 9 -l/2)- (27) 



r { LQC) = ( 1 - ^ V 1 - 3Z ) • ( 28 ) 



In the LQC model, we have the following expressions of r and s as 
and 

■Hlqc) = Z. (29) 
It is interesting to note that the pair {tlqc-, slqc} yields the ACDM (cosmological constant-cold dark matter model) 
{teg, seg} = {1, 0} when Z = or e = |. 
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Fig. 1 Fig. la Fig. lb Fig. lc 

Figs.l, la, lb, lc : The dimensionless density parameters are plotted against e-folding time. The initial 

condition is v(0) = 0.9, u(0) = 0.2. The other parameters are fixed at oj = —1,C = 1 and p\ = 10. The 
interactions are respectively b = 0.01,0.1,0.2 and 0.25 
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Fig. 2 Fig. 3 Fig. 4 

Fig 2 : The dimensionless density parameters are plotted against e-folding time. The initial condition is 

v(0) = 0.6, it(0) = 0.6. The other parameters are fixed at oj = —1, b = 0.01, C = 1 and p\ = 10. 

Fig 3 : The phase diagram of the parameters depicting an attractor solution. The initial conditions chosen 

are u(0) = 0.5, u(0) = 0.6 (green); v(0) = 0.6, u(0) = 0.6 (blue); v{0) = 0.7, u(0) = 0.6 (red); «(0) = 

0.8, u(0) = 0.6 (brown). Other parameters are fixed at to = —1,6 = 0.1, C = 1 and p\ = 10. 

Fig 4 :The deceleration parameter is plotted against the state parameter. Other parameters are fixed at 

u = -1, b = 0.01, C = 1, pi = 10, e = 0.4. 
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Fig. 5 Fig. 6 Fig. 7 

Fig 5 : The statefinder parameter r is plotted against the state parameter to gccg . The other parameters are 
fixed at oj = -1, b = 0.01, C = 1, p x = 10, e = 0.4. 

Fig 6 : The statefinder parameter s is plotted against the state parameter to gcC g. The other parameters are 
fixed at oj = -1, 6 = 0.01, C = 1, p x = 10, e = 0.4. 

Fig 7 : The ratio of density parameters is shown against e-folding time. The initial conditions chosen are 
v(0)=0.6, u(0)=0.6. The other parameters are fixed at oj = —l,b = 0.01, C = l,pi = 10. 
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4 Graphical Representation of the Phase plane Analysis 



Phase diagrams are drawn to determine the type of critical point obtained in this model. We discuss the results 
obtained in detail below: 

The dimensionless density parameters v and u are drawn in figures 1, la, lb, lc and 2. From the figures we see 
that v decreases, and u increases during evolution of the universe. This shows that the density of DM decreases while 
the density of DE increases as the universe evolves. So this result is consistent with the well known idea of an energy 
dominated universe. Figs 1, la, lb and lc are graphs depicting the same quantities with gradually increasing magni- 
tude of interaction between GCCG and DM. It is clearly evident from the figures that with the increase in interaction 
between GCCG and DM, their respective density parameters (u and v) become more and more comparable to each 
other. Thus this interacting model probably gives the best solution to the cosmic coincidence problem. Moreover 
from the figs.l, la, lb, lc and 2, it is seen that the graph of DE (green curve) does not show any urgency to creep up 
steadily in the vertical direction. This is a striking difference from the graphs obtained by Jamil et al [13] and Rudra 
et al [IJJ in case of MCG in LQC and MCG in Braneworld respectively. So in this case it is understandable that for 
GCCG, the the relative domination of the DE component over DM is far lesser as compared to MCG. The possible 
reason may be a weaker negative pressure of GCCG compared to that of MCG. This may well be the underlying- 
reason for GCCG supporting a 'No Big Rip' cosmology. The phase space diagram (figure 3) shows the attractor 
solution. The eigen values are calculated at both the critical points. At the first critical point eigen values are found 
to be (0.0381055, 1.52435). So the critical point is an unstable node. At the second critical point the eigen values are 
(0.0381047, - 6.09814). Hence the critical point is a saddle point. 

Figure 4 shows the variations of the deceleration parameter, q against w gccg . From figure 4, it is evident that 
there is a gradual decrease in the deceleration parameter g, and finally in the late universe it attains negative values, 
which suggests that there should be an acceleration in the late universe. This result is in accordance with the various 
observational data of la supernovae and CMB data which suggests that the universe is undergoing an accelerated 
expansion of late. Figure 5 and 6 shows the variations of the statefinder parameters r and s against u gcC g- From these 
figures it can be seen that r tends towards 1 and s tends towards 0. Therefore it is evident that these results tends 
towards the ACDM model. Figure 7 shows the variation of the ratio of the density parameters v and u against time. 
From the figure it is evident that the ratio of the above parameters decreases with time. So it can be concluded that 
there is a relative decrease in matter density with respect to the energy density. This is again consistent with the 
notion of an energy dominated universe. 

5 Study of Future Singularities 

It is a well known fact that any energy dominated model of the universe is destined to result in a future singularity. 
The study of dynamics of an accelerating universe in the presence of DE and DM is in fact incomplete without the 
study of these singularities, which are the ultimate fate of the universe. It is known that the universe dominated 
by phantom energy ends with a future singularity known as Big Rip [37], due to the violation of dominant energy 
condition (DEC). But other than this there are other types of singularities as well. Nojiri et al [48] studied the various 
types of singularities that can result from an phantom energy dominated universe. These possible singularities are 
characterized by the growth of energy and curvature at the time of occurrence of the singularity. It is found that near 
the singularity quantum effects becomes very dominant which may alleviate or even prevent these singularities. So it 
is extremely necessary to study these singularities and classify them accordingly so that we can search for methods 
to eliminate them. The appearance of all four types of future singularities in coupled fluid dark energy, F(R) theory, 
modified Gauss-Bonnet gravity and modified F(R) Horava-Lifshitz gravity was demonstrated in [49]. The universal 
procedure for resolving such singularities that may lead to bad phcnomcnological consequences was proposed. In Rudra 
et al [14] it has been shown that in case of Modified Chaplygin gas(MCG), both Type I and Type II singularities are 
possible. 

5.1 TYPE I Singularity (Big Rip singularity) 

If p — > oo , \p\ — > oo when a — > oo and t — > t s . Then the singularity formed is said to be the Type I singularity. 

In the present case by considering the GCCG equation of state from equation |T]) we find that there is no possibility 
for Typel singularity, i.e., Big Rip singularity. This is in absolute accordance with P. F. Gonzalez-Diaz who has 
successfully shown that by considering GCCG as the DE, Big Rip can easily be avoided, thus giving a singularity free 
late universe. 



7 



5.2 TYPE II Singularity (Sudden singularity) 

If p — > p s and p s ~ 0, then \p\ — > — oo for t — > t s and a — > a s , then the resulting singularity is called the Type II 
singularity. 

In this case we consider the equation of state for GCCG, like the previous case for our investigation. We see that 
if p — > p s and p s ~ 0, then \p\ — > for t — » t s and a — > a s . Hence there is no possibility of the type II singularity or 
the sudden singularity in case of GCCG. 



5.3 TYPE III Singularity 

For t — > t s , a — > a s , p — > oo and \p\ — > oo. Then the resulting singularity is Type III singularity. It is quite evident 
from the equation of state of GCCG that it does not support this type of singularity. 



5.4 TYPE IV Singularity 

For t — > t s , a — > a s , p — > and \p\ — > 0. Then the resulting singularity is Type IV singularity. This type of singularity 
is not supported by GCCG type DE. 

As a remark, one should stress that our consideration is totally classical. Nevertheless, it is expected that quantum 
gravity effects may play significant role near the singularity. It is clear that such effects may contribute to the singularity 
occurrence or removal too. Unfortunately, due to the absence of a complete quantum gravity theory only preliminary 
estimations may be done. 



6 Concluding Remarks 

In this work, we have considered a combination of Generalised cosmic Chaplygin gas in Loop quantum gravity model. 
Our basic idea was to study the background dynamics of GCCG in detail when it is incorporated in Loop quantum 
gravity. Dynamical system analysis had been carried out, critical points were found and the stability of the system 
around those critical points was tested. Graphical analysis was done to get an explicit picture of the outcome of the 
work. In order to find a solution for the cosmic coincidence problem, a suitable interaction between DE and DM was 
considered. Figures of density parameters were drawn for different values of interaction. It was found that increase 
in interaction resulted in more and more comparable values of the density parameters of GCCG and DM. Since the 
tendency of DE domination over DM is lesser in case of GCCG compared to MCG, GCCG is identified as a dark 
fluid with a lesser negative pressure compared to MCG. The dynamical system of equations characterizing the system 
was formed and a stable scaling solution was obtained. Hence this work can be considered to be a significant one, 
if not the best, as far as the solution of cosmic coincidence problem is concerned. Study of future singularities had 
been carried out in detail. The model was investigated for all possible types of future singularities. From the above 
analysis we conclude that the combination of GCCG in loop quantum gravity gives a perfect singularity free model 
for an expanding universe undergoing a late acceleration. 
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